Bistability in sine-Gordon: the ideal switch 
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The sine-Gordon equation, used as the representative nonlinear wave equation, presents a bistable 
behavior resulting from nonlinearity and generating hysteresis properties. We show that the process 
can be understood in a comprehensive analytical formulation and that it is a generic property of 
nonlinear systems possessing a natural band gap. The approach allows to discover that sine-Gordon 
can work as an ideal switch by reaching a transmissive regime with vanishing driving amplitude. 
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I. INTRODUCTION 

A nonlinear medium submitted to wave irradiation at 
a frequency in a forbidden band gap can undergo bistable 
behavior and present hysteresis properties. This bistabil- 
ity has attracted much attention, e.g. in nonlinear optics 
as a means for a medium to switch from total reflection 
to partial (sometimes total) transmission []J, or in su- 
perconducting junction devices as a means to conceive 
amplifiers that "remain efficient in the quantum limit" 
0. 

We attempt to give here a comprehensive interpreta- 
tion of this phenomenon in terms of both analytical de- 
scription and numerical simulations, in order to unveil 
a particular stationary regime presenting a non-zero out- 
put for vanishing input, what we call the ideal switch and 
which allows for detection of (almost) vanishing signal. 

To that end we consider the sine-Gordon equation on 
the finite interval x S [0, L] 



u tt 



0, 



associated to the boundary value problem 

u(0,t) = /(*), «xCM) = o, 



(1) 



(2) 







on a vanishing initial state (with /(0) = and /t(Q) : 
for compatibility of the initial and boundary values). 

This is a quite standard problem in physics of a wave 
equation with a forced extremity (Dirichlet boundary 
condition in x — 0) and a free end (Neumann boundary 
condition in x = L), associated to Cauchy initial data at 
t = 0. A related physical situation is for instance a long 
Josephson junction |3, Q or an array of coupled short 
junctions (Josephson superlattice) 0,0. Note that, de- 
pending on the used external driving, the boundary @ 
has possibly to be replaced with u x (0, t). 

An important subclass of boundary f(t) is constant 
amplitude periodic driving at a frequency in the natural 
band gap of the system, namely 



f(t) = Bo cos(fit), ft < 1, 



(3) 



tial shock. While for a linear system this boundary exci- 
tation does not flow through, nonlinearity allows for en- 
ergy transmission above the threshold amplitude which 
reads (in the semi-infinite case L — ► oo) 



B s = 4 arctan(6 s 



b = 



l-O 2 



(4) 



after a convenient transient sequence where f(t) grows 
from a vanishing amplitude to the value Bq to avoid ini- 



This is called nonlinear supratransmission [7] and hap- 
pens by emission of solitons (moving breathers) that 
propagate in the nonlinear medium. 

This process, quite generic, has been experimentally 
realized on a chain of coupled pendula and applies 
for instance in discrete systems of coupled waveguide ar- 
rays [jj where the forbidden gap results from discreteness, 
or else in Bragg media (periodic dielectric structures) un- 
der constant micro- wave irradiation in the photonic band 
gap [T3- In Josephson junctions arrays, submitted to mi- 
crowave excitation, the boundary u x (0, t) = f(t) induces 
the threshold B s = 2(1 - tt 2 ) [J. 

In the finite line case considered here, we shall again 
find a nonlinear supratransmission threshold which tends 
to the value J3J f° r large L. But a property far less un- 
derstood is the hysteresis loop obtained by decreasing 
the amplitude excitation Bq from the threshold B s . This 
property has been for instance observed on numerical 
simulations 0, |(| in the context of Josephson superlat- 
tices, but both the analytical expression of the threshold 
and the very nonlinear mechanism involved have not been 
clarified. 

We shall establish a general procedure to determine 
the threshold by studying the standing periodic solu- 
tions of sine-Gordon which synchronize to the driving 
frequency il and adapts to the driving amplitude Bq. Al- 
though these two conditions are sufficient to determine 
completely the solution, it is not uniquely defined. In- 
deed we shall prove that a fixed set of physical parame- 
ters {L, Q, i?o} may be related to more than one solu- 
tion. This is the principle that leads to bistability when 
B < B s . 

As an interesting consequence we obtain that there ex- 
ists a regime where a vanishing input amplitude Bq pro- 
duces a non-vanishing output amplitude. This process 
shows that sine- Gordon can be thought of as an ideal 



2 



switch along the hysteresis loop from zero to zero input 
amplitudes. 

The paper is organized as follows: in the next section 
we display the set of explicit solutions to sine-Gordon on 
a length L submitted to the only requirements that the 
input boundary amplitude be Bo and the period of the 
solution be 2w/fl. The following section is devoted to 
the analytical definition and evaluation of the threshold 
of bistability. Then we show by numerical simulations in 
section 4 that those explicit solution are indeed produced 
by the boundary driving @ and we check bistability pre- 
dictions. In particular we compute, for a more realistic 
damped sine-Gordon model, the power released by the 
boundary driver to the medium and find that after hav- 
ing switched, this power is 2 to 3 orders of magnitude 
larger than before the switch. 



The second fundamental hypothesis consists in express- 
ing that the solution adapts to the driving amplitude 
Bq = 4 arctan(a) , which gives 

X(0) = a. (12) 

The two relations (| 1 1 11 and (|12|l constitute a closed sys- 
tem of equations for the two unknowns A and a in terms 
of the physical constants a, f2 and L. 

The point is that bistability occurs because the solu- 
tion of ©© drastically depends on the sign of a. We 
shall indeed discover that there may exist different solu- 
tions that do synchronize to f2 and adapts to a. In other 
words, for any fixed Q and L, a given input amplitude 
a may correspond to more than one value of the output 
amplitude A as depicted on figH 



II. EXPLICIT SOLUTIONS 



A. General expressions. 

Under boundary condition (|3J) , in order to describe the 
periodic asymptotic regime reached in numerical simula- 
tions, we follow [ll| an d seek a solution 

u(x,t) = 4arctan[A(x)T(t)]. (5) 

The boundary condition in u x {L,t) — then reads 

X'(L) = 0, X(L) = A, (6) 

where we have defined the amplitude parameter A such 
as to scale T{t) to unity, in other words 



3t : T'(i ) = 0, T(t ) = l. 



(7) 



By inserting expression (JSJ) in the sine-Gordon equation 
Q and by use of constraints © and Q , we obtain dif- 
ferential equations with a unique free parameter a: 

(X') 2 =aT{A 2 -X 2 ){X 2 + ^), (8) 



(T') 2 = a(l - T 2 )(T 2 + T), 
where prime denotes differentiation and where 



r = 4 



l 



A 2 a(l + A 2 )' 



(9) 



(10) 



Thanks to (jfJJl and Q) the equation for X(x) is integrated 
on [x, L] and the one for T(t) on [t, t ] ■ The solution is 
then completely defined (in terms of elliptic integrals) by 
the values of the two parameters A and a, determined as 
follows. 

Our first fundamental hypothesis is to assume, accord- 
ingly with numerical simulations, that the solution syn- 
chronizes to the boundary driving, namely that the func- 
tion T(t) is periodic with the period of the driver: 



T(t- 



2tt, 

n 1 



T(t). 



(11) 



B. Type I solutions. 

We call type I solutions those obtained for a > (im- 
plying r > 0) for which we obtain [l2l | 



T(t) = cn(w(t-t ),i/), 
X(x) = A cn(k(x - L),/x), 

uj 2 =a(l + T), v 2 = 



1 + V 



TA 2 ' 



1 + TA 4 ' 



(13) 
(14) 



(15) 



where cn(-,m) is the cosine- amplitude Jacobi elliptic 
function of modulus m. According to the result- 
ing solution u(x, t) is called plasma oscillation and we 
have from l|15|l 



CO 



I- A 2 
l + A 2 ' 



(16) 



This relation between the nonlinear wave parameters to 
and k is often called a nonlinear dispersion relation but 
we shall reserve such a denomination to the true disper- 
sion relation which relates the actual period 4K(fi) /k of 
X{x) to the period 4K(v)/io of T(t). 

The parameters a and A are now determined by requir- 
ing synchronization (|llfl and input datum (|12|) . namely 
here by solving for unknowns a and A the system 



n k(i>) = - u, 



A cn(kL, //). 



(17) 



(Note: the complete elliptic integral K(^) is well defined 
as for T > we have < v 2 < 1.) 

It is useful to express the above solution in terms of 
the two wave parameters lu and k. This is done by using 
(|16fl to eliminate A and a from the definitions of /i and 
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v. We get 



1-uj 2 



v 2 = 



l+UJ 2 -k 2 ' 

k A - {uj 2 - I) 2 



4w 2 
4fc 2 



(18) 



(1 + fc 2 ) 2 -w 4 



System (11711 appears then as an equation for the deter- 
mination of the parameters uj and k from the data of the 
length L, the boundary driver's frequency ft and ampli- 
tude a. 



C. Type II solutions. 

New types of solutions are obtained for a < for which 
the evolution (0 of T{t) requires T < in order to guar- 
antee the constraint Q). Defining then 







-a. 



A 



1 



13(1 + A 2 ) A 21 
the basic equations (JSJ and © become 
(T 1 ) 2 = [3{l-T 2 )(K-T 2 ), 
(X') 2 = ?K(A 2 - X 2 )(X 2 - 



(19) 

(20) 
(21) 



It appears that the constraint Q which states that 
T(io) = 1, requires A > 1, namely 



0< 



A 2 



(1 + A 2 ) 2 ' 



(22) 



a condition that must be checked a posteriori when com- 
puting j3 from the synchronization constraint. 

The structure of the equation (|2~TJ) implies two classes 
of solutions depending on the relative values of A 2 and 
1/(AA 2 ). Type II solutions are obtained for A A 4 > 1 
which, together with constraint l|22|l . reads 



A 2 > 1 
A 2 < 1 



< [3 < A 2 /(l + A 2 ) 2 , 
0</?< A 4 /(l + A 2 ) 2 , 



The solution of (|20|H|21II can now be obtained as 

T(t) = sn(w(t -ti),z/), ti = t + K(i/)/w 
X(ai) = A dn(k(x-L),/i), 
1 



" = A' 



w 2 = /3A, 
fc 2 = /?AA 2 , ^ 2 = 1 
and we have the relation 

J 1 



AA 4: 



A 2 ' 



(23) 
(24) 



(25) 
(26) 



(27) 



(28) 



between the wave parameters k and u> for type II solution. 

The parameters /3 and A are determined as before 
by requiring synchronization (fTTJl and input datum ltl"2*|) . 
namely 



ft K(t/) = - uj, 



A dn(kL,/i). 



(29) 



Although the above equation, for real valued parameters, 
has two sets of solutions {(3, A}, only one set verifies the 
constraint A A 4 > 1. 

As before, we express the type II solution in terms 
of the wave parameters u> and k, by means of 128|) to 
eliminate A and j3 in fi and v. We obtain 



k 2 

A 2 - — 

2 ' 

2 k 2 1 

7/ Z = 

,2 



■ UJ 



k 2 



uj 2 uj 2 + k 2 



(30) 



^ = 1- 



1 - uj 2 - k 2 



k 2 



k 2 



System (12911 determines then the parameters uj and k 
from the data of L, ft and a. 
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FIG. 1: Plot of the curves |A(a)| in the three cases for L — 
5 and Q — 0.6. The crosses indicate the points where the 
solution changes from one type to the other. The vertical line 
shows the threshold amplitude a s (next section) above which 
supratransmission occurs by emission of solitons. 



D. Type III solutions. 

The type III solution is obtained still for a < when 
A A 4 < 1. Such can be realized only in the case A 2 < 1 
by requiring 



A 4 /(l + A 2 ) 2 <p< A 2 /(I + A 2 ) 2 . 

The solution of J3TJ now reads 

T(t) = sn(w(t - ti), v), ti = t + K(v)/u 
X(x) = Adn _1 (k(x-L),^), 

u 2 = PA, ^ 2 = {, k 2 = JL, M 2 = 1-AA 4 



(31) 

(32) 
(33) 

(34) 
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and the wave parameters obey 

1 



uj 



1 + A 2 



(35) 



The synchronization condition and input datum 1121) 
furnish here the system 

ft KM = £ uj, a = —4 r- (36) 

for the unknowns parameters (3 and A. The same remark 
as for the type II solution holds here, namely that the 
system (|36|l has two sets of solutions ^4} but only one 
verifies the constraint KA A < 1. 

Again the type III solution can be expressed entirely 
in terms of the two parameters u> and k by means of l|35|l 
which gives 



A l = 



v 2 = 



l-LU 2 -k 2 

uj 2 + k 2 



j 2 + k 2 ' 



(37) 



M 1 — To 



1-LO 2 - k 2 



Then l|36ll is a system for the parameters uj and k in terms 
of the data L, fl and a. 



III. BISTABILITY THRESHOLDS 

The above 3 solutions are now used to describe analyt- 
ically bistability properties of sine-Gordon. The first step 
is to define and calculate, at given length L, the thresh- 
old a s as functions of the driving frequency f2. Last, 
decreasing the input amplitude from a s , once a trans- 
mission regime has been reached, the system locks to the 
type-I solution which holds down to a vanishing driving 
amplitude. This is a property that makes sine-Gordon as 
the ideal switch and allows to understand how it can be 
used to amplify weak (vanishing) signals. 

A. Transmission threshold 

As shown by figure ^ increasing the input amplitude 
from a = generates the type III solution. This solu- 
tion has a maximum input value a s resulting from l|36|) 
as the point where dn(kL, /i) reaches its minimum value 



y/l — /i 2 , namely 



a 2 = 



.4 2 



(38) 



Such a definition of the threshold a s is more conveniently 
written in terms of uj and k through l|37|) as 



*1 



(39) 



In this equation, the parameters uj and k are determined 
through l|36|) . which, at the threshold, can also be written 
as the system 



QK(u) = -oj, kL = K(p), 



(40) 



with v and /i given by H37J) . 




0.2 0.4 0.6 0.8 

DRIVING FREQUENCY 



FIG. 2: Input amplitude 4arctan(a s ) at the threshold, solu- 
tion of 1391 14UH as a function of SI for L = 5 (crosses) and 
compared to its limit value 4arctan(6 s ) as L — > oo, given by 
<H3ll (full line). 

The figure shows that the threshold a s is quite close 
to the expression of b s in Q). This property is demon- 
strated in general by studying the limit L — » oo in the 
type III solution. At threshold amplitude, the relation 
(I40|l between L and /x gives the necessary condition 

L^oo fj,-yl cj 2 ^l-fc 2 , (41) 

with which the synchronization condition provides 

L^oo => v ^ => uj ^Q. (42) 

With this in hands, the expression Q39J1 of the threshold 
readily gives the expression namely 

1-9? 



L 



O 2 



(43) 



Remark: it is instructive to compute also the limit 
L — > oo on the solution itself at the threshold where 
from l|37|l and (|41|l obviously A — > 0. In that case we 
rewrite the solution X{x) of ll3*3l) as 

a dn(kL, [i) 



X(x) 



(44) 



dn(k(x — L), /i) ' 

expand the denominator, take the limit L — > oo first and 
make then fj, —> 1. We obtain X(x) of (|33|l as 

X(x) — > a s sech(kx). (45) 



The same procedure applied to T(t) in l|32() . with ^ 
provides 



0, 



r(t) 



cos(w(t - 1 )), 



(46) 



and the resulting solution u(x,t) of sine-Gordon on the 
semi-line is the stationnary breather centered in x = 
accordingly with ■ 
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IDEAL SWITCH (L^5) 




0.2 0.4 0.6 0.8 
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FIG. 3: Output amplitude 4arctan(^4) when a = as a func- 
tion of SI for L — 5 obtained from (I48H by solving 1471 . 



We observe that, for a given length, there exists a 
threshold in frequency below which no ideal switching 
is allowed. This is understood by observing that A 2 di- 
verges when uj 2 — > k 2 — 1, for which fj, — * 1 (a limit 
threshold that has exactly the same origin as in (|41(l when 
L —> oo). Conversely, at given driver frequency Q, there 
exists a minimum length L of the medium to obtain an 
ideal switch, it is displayed on fig. ^ 

Let us remark that the notion of nonlinear dispersion 
relation is not useful to predict, at given driving fre- 
quency f2, the minimum driver amplitude that generates 
transmission, as indeed we have here an example where 
this minimum is simply vanishing. 



IV. NUMERICAL SIMULATIONS 



B. Ideal switch 



A. Damping and boundary driving. 



Considering the type I solution, expression (|17|l that 
links the input a to the output A can produce a = with 
A / 0, such as to generate a regime of non- vanishing 
output value with a vanishing input amplitude, the ideal 
switching regime. This is the case when 



Bistable properties, and in particular ideal switching, 
have been analytically described in the integrable case 
(Q. However, any realistic physical situation must take 
into account the damping inherent to the medium. The 
simplest way to include damping is to assume the model 



kL = K(fi), n K(u) = - lu. 



(47) 



where A, /i and v are related to w and k through (|18|l . 
Note the formal analogy with l|4U|) where the parameters 
A, fi and v are different. 




0.2 0.4 0.6 

DRIVING FREQUENCY I 



smu = 0, 



(49) 



associated with the initial-boundary value problem J5J in 
the particular subclass J3J. We study here the bistable 
properties of l|49(l by numerical simulations and compare 
the results to the analytical predictions. 




FIG. 4: Minimal length of the system below which one does 
not have ideal switch (the solution of l|47^ ceases to exist). 

This is a system of equations for {ui, k} whose solution 
then produces the seeked output amplitude by 



A 2 



1 



1 + UJ 2 - k 2 ' 



(48) 



defined by (|47|l in terms of the physical entries L and £1. 
We have plotted in figOHthe output amplitude 4arctanA 
in the ideal switching case (a = 0) as a function of $1 for 
length L = 5. 



FIG. 5: Two different paths for driving the sine-Gordon sys- 
tem. Bo is the driving amplitude in stationary regime and B s 
is a supratransmission threshold. 

To that end the system is driven at the input boundary 
x = with a bandgap frequency and time dependent am- 
plitude as follows: in a first numerical simulation, the am- 
plitude is smoothly increased up to the value Bq smaller 
than the supratransmission threshold B s . In a second 
simulation, the amplitude is increased up to a value ex- 
ceeding the supratransmission threshold and then, after 
a time sufficient to generate moving breathers, it is de- 
creased to the same value Bq as in the first case. The 



TIME [t] g50 1 DISTANCE [X] 




950 

FIG. 6: Three dimensional plots of the sine-Gordon system 
dynamics after having reached stationary regime. Upper plot 
corresponds to the driving path 1 of fig|^]and lower plot to 
path 2. In both cases driving frequency is Q = 0.5, stationary 
driving amplitude Bo = 0.3, damping coefficient 7 = 0.01 and 
system length L — 4.1. 

figure |S] displays the two time variations of the driving 
amplitude that we have used in the numerical simula- 
tions. After having reached a stationary regime, although 
the driving amplitudes B are equal in both cases, the 
dynamics drastically differ in those two cases as it is ex- 
pected from the analytical consideration presented above 
and described hereafter. 



FIG. 7: Comparison of analytical expressions (solid lines) for 
the profiles of standing waves 1261 . CIJ, corresponding 

to solutions of type III, II and I respectively, and numerical 
simulations (circles). The type III solution is reached along 
path I of fig. Awhile the type I results from path 2. The type 
II solution, unstable, is never reproduced by the system. 



system has locked to a stationary solution with a small 
driving amplitude (here Bq = 0.3) and a large output am- 
plitude (evaluated at B = 4.35). Let us mention that we 
can drive the system with amplitudes down to Bq = 0.1 
and still have the type I solution (large output ampli- 
tude) despite presence of damping in the system, getting 
thus a regime of almost ideal switch, or almost perfect 
detector. 

Another remark is that the system never locks to the 
exact solution of type II (|25[lil2rj[) . simply because this 
solution is unstable. For instance, if that exact solution is 
used as an initial condition in sine-Gordon, it eventually 
breaks down. 

To be complete, we compare the analytical expressions 
of the standing waves profiles of l|3lfl) . l|2^|) and iTr4"|) to 
the results of numerical simulations in the case of the two 
different driving paths in the context of figEI The result 
is displayed in figure|Z|where the obtained perfect match- 
ing shows that indeed the system locks to the analytical 
solution obtained by assuming frequency synchronization 



B. Evidence of bistability. 

In most of our numerical simulations we choose the 
driving frequency in the middle of the band gap ft = 
0.5, use damping parameter 7 = 0.01 and length L = 
4.1. For the driving amplitude along path 1 in figEJ we 
always observe a stationary regime with decaying profile 
of the standing wave, very well described by the exact 
analytical solution of type III l|32|ll|33() . Instead, when 
driving along path 2, we get the picture corresponding to 
the exact type I solution l|13 |) l|14 |) . Those two drastically 
different behaviors of the system are displayed as three 
dimensional plots in figE) 

It is remarkable on the second picture of figElthat the 




INPUT AMPLITUDE 

FIG. 8: Comparison of input-output amplitude dependence 
obtained from numerical simulations (circles) with analytical 
curves derived from formulas 1171 . 12911 . 1361 . 



7 




1 2 3 

DRIVING AMPLITUDE 



FIG. 9: Power P of f r >(H in terms of the driving amplitude in 
lin-log (main plot) and lin-lin (inset) scales. Solid line (solu- 
tion I|13ll4fl ') and circles (simulations) correspond to driving 
regime 1 of figEl Dashed line (solution 132I33H ) and asterisks 
(simulations) result from driving regime 2. Length is L = 4.1, 
driving frequency SI = 0.5 and damping coefficient 7 = 0.01. 



and amplitude matching. This is completed by plotting 
in fig|S]the input-output amplitude dependence obtained 
from numerical simulations and its comparison with an- 
alytical curves derived from formulas l|17fl . (|29fl and (1301 . 



C. Energetic considerations. 

The physically useful bistable nature of the system 
manifests in large difference between the energy dissipa- 
tion in the two stationary regimes. The averaged energy 
released from the driver in unit time, can be expressed 
as 



n 



2tt/0 L 



P = 2^ J dt J dxiv %- 




(50) 



Indeed one can easily obtain 

L L 



J dxH = - J dx~fu 2 t - {u x u t )\ x=0 , 



1 , 1 2 

n = —u t H — uz + 1 — cos it, 

2 2 



(51) 



which by averaging on one period furnishes (|50() as the 
boundary value in x = is periodic. 

Although the analytical solution of the sine-Gordon 
equation are not solutions of the damped version 149|l , 
we may compare the power defined above obtained by nu- 
merical simulations of H49|) to the expression Q5l)f) where 
u is simply replaced by the exact solution (type III before 
the switch, type I after). The result of this comparison is 
displayed in figElwher we see that expression (|50|) with 
analytical solutions fit strikingly well the numerical sim- 
ulations, and that the power P after the switch is two to 
three orders of magnitude greater than before the switch. 
V. COMMENTS AND CONCLUSION. 



It is worth mentioning that while the analytical solu- 
tion of type I holds for any length L, in a realistic physi- 
cal system (nonzero damping), the situation is different. 
In particular, for large L and low driving amplitudes, 
when several nodes of type I standing wave solution are 
present, the solution cannot survive and decays to type 
III solution. This is understood by the following sim- 
ple argument: all of the exact solutions derived in the 
previous sections are standing waves, i.e. they do not 
generate energy flux. Thus, regions of the system far 
from the boundary cannot gain energy from the driver 
and the oscillations will eventually fade away. 

We have essentially demonstrated, both by analytical 
and numerical treatment, that the bistable property on 
the sine-Gordon system allows to generate a particular 
regime that works as an ideal switch: nonzero output for 
vanishing input. In a realistic physical system (including 
damping) the property is conserved but for small (non- 
vanishing) input. This nonlinear hysteretis has been 
shown to correspond to quite determinant differences in 
the power released by the driver to the system. 
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